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Abstract 



Adaptive dynamics is a widely used framework for modeling long-term evolu- 
tion of continuous phenotypes. It is based on invasion fitness functions, which de- 
termine selection gradients and the canonical equation of adaptive dynamics. Even 
though the derivation of the adaptive dynamics from a given invasion fitness func- 
tion is general and model-independent, the derivation of the invasion fitness function 
itself requires specification of an underlying ecological model. Therefore, evolution- 
ary insights gained from adaptive dynamics models are generally model-dependent. 
Logistic models for symmetric, frequency-dependent competition are widely used 
in this context. Such models have the property that the selection gradients derived 
from them are gradients of scalar functions, which reflects a certain gradient property 
of the corresponding invasion fitness function. We show that any adaptive dynam- 
ics model that is based on an invasion fitness functions with this gradient property 
can be transformed into a generalized symmetric competition model. This provides 
a precise delineation of the generality of results derived from competition models. 
Roughly speaking, to understand the adaptive dynamics of the class of models sat- 
isfying a certain gradient condition, one only needs a complete understanding of the 
adaptive dynamics of symmetric, frequency-dependent competition. We show how 
this result can be applied to number of basic issues in evolutionary theory. 



2 



1 Introduction 



Adaptive dynamics ( |Metz et aL\ ( |1996[ ), |Geritz et al. \ ( |1998| ), pieckmann & Law 



(1996)) has emerged as a widely used framework for modeling long-term evolution 
of continuous phenotypes. The basic ingredient of an adaptive dynamics model is the 



invasion fitness function (Metz et al. ( 1992)), which describes the ecological growth 
rate of rare mutant phenotypes in a given resident community, which is assumed to 
persist on a community-dynamical attractor. The invasion fitness function determines 
the selection gradients, which are in turn the core ingredient for deriving the canoni- 
cal equation ( pieckmann & Law ( 1996 )) for the adaptive dynamics of the phenotypes 
under consideration. Following this basic recipe, adaptive dynamics models have 
been constructed for a plethora of different ecological settings, and have been used 
to analyze a number of interesting and fundamental evolutionary scenarios, such as 



evolutionary cycling in predator-prey arms races (Dieckmann et al. ( 1995 ), Marrow 



et al.\ ( [199 6)), evolutionary diversification (Dieckma nn & Doebeli| ([1999), Dieck- 



mann et a7] ( |2004[ ), poeberi| ( |201 1| )) and evolutionary suicide Gyllenber g & Parvinen 



(2001 1, Parvinen (2005 ). Even though the derivation of the adaptive dynamics from 
a given invasion fitness function is general and model-independent, the derivation of 
the invasion fitness function itself requires specification of an underlying ecological 
model. Therefore, evolutionary insights gained from adaptive dynamics models are 
generally tied to a specific ecological setting, and hence model-dependent. 

One particular ecological model that has been often used to derive invasion fitness 
functions and adaptive dynamics is the symmetric logistic competition model, which 
in fact is the most popular ecological model in the theory of ecology and evolution. 
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The basic form of this model is 

dN ( N\ 

rN(l--), (1) 



dt V K . 

where N is population density, and r and K are parameters describing the intrinsic 
per capita growth rate and the equilibrium population size, respectively. K is often 
called the carrying capacity of the population, but it useful to note that K can also 
be interpreted as a property of individuals, i.e., as a measure of how well individuals 
cope with competition, as expressed in the per capita death rate rN/K. The eco- 
logical model (1) can be used to construct a well-known adaptive dynamics model 
by assuming that the carrying capacity K(x) is a positive function that depends on 
a continuous, 1-dimensional phenotype x (e.g., body size), and that the competitive 
impact between individuals of phenotypes x and y is given by a competition kernel 
a(x,y). For simplicity, it is often assumed that the intrinsic growth rate r is inde- 
pendent of the phenotype x, and hence is set to r = 1. To derive the corresponding 
invasion fitness function, it is assumed that a resident type x is at its ecological equi- 
librium density K(x). The dynamics of the density N(y) of a mutant type y is then 
given by 

dN(y) r _ a(x,y)K(x) + a(y,y)N(y) \ (2) 



dt vyy V K(y) ) 

where a(x, y)K(x) is the competitive impact that the resident population exerts on 
mutant individuals. Assuming that the mutant is rare, N(y) « 0, the per capita 
growth rate of the mutant type y in the resident x is 

/(,,,) = l- gfejgW (3) 



This is the invasion fitness function for the given ecological scenario (Doebeli (201 1 )) 
The corresponding selection gradient is 

df(x,y) 



dy 



(4) 

y=x 



The selection gradient in turn is the main determinant of the adaptive dynamics. 
More precisely, for 1 -dimensional traits x the adaptive dynamics is 

^ = m{x) ■ s(x), (5) 

where m(x) > is a scalar quantity that reflects the rate at which mutations occur 
dMetz et Q/.| ( [T996l ), [Geritz et a/.| ( fT998| ), |Dieckmann & Law| ( |T996| )). 



In particular, singular points, i.e., equilibrium points of the adaptive dynamics, 
are solutions x* of s(x*) = 0. Symmetric competition models are characterized by 
the assumption that da(x,y)/dy\ y=x = for all x, it is also usually assumed that 
a(x, x) = 1. In this case, s(x) = K'(x)/K(x), and singular points x* are given as 
solutions of K'(x*) = 0, and hence are given as local maxima and minima of the 
carrying capacity. Moreover, singular points that are local maxima are attractors of 
the adaptive dynamics, i.e., ds/dx(x*) < 0. For logistic models it is often assumed 
that K(x) is unimodal, attaining a unique maximum at some trait value x . K(x) 
then represents a (global) stabilizing component of selection for x , and the adaptive 
dynamics (5) converges to x . 

In contrast, the competition kernel a(x,y) generally represents the frequency- 
dependent component of selection. For symmetric competition it is usually assumed 
that the effect of competition decreases with increasing phenotypic distance \x — 
y\, which implies negative frequency dependence, because it confers a competitive 
advantage to rare phenotypes. However, the opposite is also possible, so that a(x, y) 
has a local minimum as a function of y at y = x, which can occur for example in 



models with explicit resource dynamics (Ackermann & Doebeli (2004)). While the 



carrying capacity is the sole determinant of the singular points and their convergence 
stability for the adaptive dynamics (5), the competition kernel comes into play when 
determining the evolutionary stability of the singular point x*. Evolutionary stability 
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of x* is determined by the second derivative of the invasion fitness function at the 
singular point: 



d 2 f(x%y) 



dy 2 



y=x* 



K"{x*) _ d 2 a(x*,y) 
K(x*) 



dy 2 



(6) 



y=x' 



The singular point x* is evolutionarily stable if and only if expression (6) is negative, 
and it is clear that if d 2 a(x*, y)/dy 2 \ y=x * is negative enough, then this condition will 
not be satisfied, and instead x* will be evolutionarily unstable. In particular, the dis- 
tinction between convergence stability and evolutionary stability makes it clear that 
it is possible for the singular point x* to be both convergent stable and evolutionarily 
unstable. In this case, x* is called an evolutionary branching point, for such points 
are potential starting points for evolutionary diversification. The phenomenon of evo- 
lutionary branching is an iconic feature of adaptive dynamics and has been studied 



extensively ( |Geritz et a/.] ( [T9^ , |Doebeli| ( |20lT] )). 

It is straightforward to extend the symmetric logistic competition model to multi- 
dimensional phenotype spaces ( |Dieckmann & Law] ( |1996[ ), |Leimar| ( |2009| ), |Doebeli 
& Ispolatov (2010), Doebeli ( 201 ![ )). In this case, x E W n is a m-dimensional 
vector, where m is the dimension of phenotype space, and K(x) : IR m — > R is a 
scalar function, as is the competition kernel a(x, y) : W n x IR m — >■ BL For symmetric 
competition, it is assumed that the partial derivatives da(x,y)/dyi\ y=x = for all 
x and all i = 1, ...,m (as well as, without loss of generality, ) = 1 for all 

x). The corresponding invasion fitness function f(x, y) : M m x M. m — > M. again has 
the form (3), and the selection gradient is given as a vector- valued function s(x) = 

OiO), . . . , s m (x)), : K m R, where 
df(x,y) 



SAX 



dyi 



i — 1, m 



(7) 



y=x 



The adaptive dynamics is then given by a system of m coupled differential equations 

dx 
~dt 



M(x) ■ s(x) 



(8) 



where M(x) is the mxm mutational variance-covariance matrix, specifying the rate 
and magnitude at which mutations occur in the various trait components Xi, as well 



how mutations in different trait components are correlated (Leimar (2009), Doebeli 



(201 1 )). The matrix M(x) is typically assumed to be symmetric and positive definite. 

For symmetric competition models, it is easy to see that the m components of the 
selection gradient are 

1 dK{x) 



Sl (x)= df{x > y) 



K(x) dxi 



(9) 



for i — 1, m. Thus, the selection gradient is the gradient of a scalar function, 

s(i) = VS(i), (10) 

with S(x) = \n(K(x)). Thus, just as in the 1-dimensional symmetric competi- 
tion model, the selection gradients in multi-dimensional generalizations of symmet- 
ric competition are essentially gradients of the carrying capacity function K(x) (i.e., 
gradients of the stabilizing component of selection). The selection gradients there- 
fore induce an evolutionary hill-climbing process towards local maxima of the car- 
rying capacity. The adaptive dynamics (8), resulting from applying the mutational 
variance-covariance matrix to the selection gradient, is a "warped" version of the 
hill-climbing process generated by the selection gradients. If the mutational matrix 
M{x) is positive definite, this warped hill-climbing process is essentially equivalent 
to the unwarped version defined by the selection gradients alone. In particular, in this 
case the adaptive dynamics (8) also converges to local maxima of K(x). 

It is worth noting that the structure of the canonical equation (8) is similar to other 



general equations for evolutionary dynamics, such as those introduced by Lande 



(1979) and those introduced by |Nowak & Sigmund (1990). But the assumptions 



underlying the canonical equation (8), and in particular the notion of invasion fit- 
ness, are unique features of adaptive dynamics. In general, invasion fitness functions 



f(x, y) can be derived for multidimensional phenotypes in each of a number of in- 



teracting and coexisting species for many different ecological scenarios (Dieckmann 



& Law ( 1996), Leimar (2009), Doebeli (2011 )). In each of the interacting species, 



invasion fitness is the long-term per capita growth rate of rare mutant types y in a 
population in equilibrium that is monomorphic for the resident type x. Here we con- 
sider the following question: under what conditions does an arbitrary invasion fitness 
function for multi-dimensional phenotypes in a single species has the form of an in- 
vasion fitness function derived from a logistic symmetric competition model? This 
is a relevant question because symmetric competition models have been used for a 
long time as a basic metaphor to generate ecological and evolutionary insights. It is 
therefore of interest to understand how universal such models are. For example, in 



Doebeli & Ispolatov (2010) it has recently been shown that evolutionary branching, 



and hence adaptive diversification, becomes more likely in symmetric competition 
models if the dimension of phenotype space is increased, and it would be useful to 
know whether this applies to other models. In the present paper we give a general 
and precise condition for any given invasion fitness to be equivalent to the invasion 
fitness function derived from a symmetric competition model. We also show how this 
result can be applied to shed light on some general issues in adaptive dynamics the- 
ory, such as the notion of frequency-dependent selection, the relationship between 
symmetric and asymmetric competition, the existence of complicated evolutionary 
dynamics, and the problem of evolutionary stability in single- species models. 
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2 A condition for universality 



We consider an invasion fitness function f(x,y) for a single species, in which the 
multidimensional trait x = (xi, . . . ,x m ) denotes the resident trait, and the vector 
y = (yi, . . . , y m ) denotes the mutant trait. Recall that the selection gradient s(x) is a 
vector 

s(x) = {si(x), s m (x)), (11) 

where 

Si(x) = B J^A (12) 



y=x 



dyi 

Dieckm ann & Law| ( |1996[ ), [Geritz et q/.] ( |1998[ ). Thus, for a given resident trait vector 



x, the selection gradient s(x) is the gradient of the fitness landscape given by f(x, y), 
but it is important to note that s(x), which is a vector- valued function on an m- 
dimensional space, is defined as the gradient with respect to y of a scalar function 
defined on a 2m-dimensional (x, y) space. In particular, the vector field s(x) is in 
general not the gradient field of a scalar function defined on m variables. 

The latter case is captured in the following Definition: We call the selection gra- 
dient s(x) a gradient field if it can be obtained as the derivative of a scalar function 
S(x), i.e., if there is a function S(x) such that 

s(x)=VS(x), (13) 

i.e., such that 

dS(x) 

Si(x) = -T7 — for i = l,...,m. (14) 

OXi 

The following proposition describes the conditions under which a given invasion 
fitness function is equivalent to an invasion fitness function derived from a symmetric 
competition model. 



Proposition 2.1 The selection gradient s(x) of an invasion fitness function f(x,y) 
is a gradient field if and only if the invasion fitness function is of the form 

/(l , v ) = i-«^M (15) 

K{y) 

for some scalar functions K(x) anda(x, y), such that K(x) > 0, a(x, x) = I for all 
x, and da(x,y)/dy\ y=x = for all x. 

Proof If the invasion fitness function f(x, y) has the form (15), let S(x) = In K(x). 
Then it is easy to calculate that 

six) = = VS(x), (16) 

K{x) 

where V is short for (d/dxi, . . . , d/dx m ). Hence the selection gradient is a gradient 
field. 

Conversely, if the selection gradient is a gradient field, s(x) = VS(x), then let 
K{x) = exp [S(x)\ , (17) 

and 

a(x,y) = [l-f(x,y)}^\. (18) 

Note that K(x) > and 6t{x,x) = 1 for all x, because f(x,x) = for any inva- 
sion fitness function (i.e., the long-term per capita growth rate of individuals with the 
resident phenotype must be 0, because the resident, assumed to exist on an equilib- 
rium, neither goes extinct nor increases without bounds). Moreover, using the fact 
that s(x) = VS(x) = VK(x) ■ K(x)~ 1 by construction, it is easy to check that 
da(x, y)/dy\ y=x = for all x, and that expression (15) holds for the invasion fitness 
function f(x, y). | 
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Thus, any invasion fitness function for which the selection gradient is a gradient 
field can be viewed as the invasion fitness function of a generalized Lotka-Volterra 
model for symmetric, frequency-dependent competition. As mentioned in the previ- 
ous section, in this case the adaptive dynamics is a hill-climbing process on a fixed 
fitness landscape determined by the (generalized) carrying capacity K(x) defined 
by (17), and singular points of the adaptive dynamics are given as local extrema of 
K(x). In particular, a singular point of the adaptive dynamics is convergent stable if 
and only if it is a local maximum of K(x). 

Next, we derive a general condition for a selection gradients to be a gradient 
field. For this we recall that a geometric object is called simply connected if any 
closed path can be shrunk continuously to a single point within that object. Thus, 
roughly speaking an object is simply connected if it doesn't have any holes. A typical 
example is Euclidean space R n for any n. 



Proposition 2.2 Assume that the phenotype space (i.e., the set of all attainable, m- 
dimensional phenotype vectors x) is simply connected. Then the selection gradient 
s(x) of an invasion fitness function f(x, y) is a gradient field if and only if the inva- 
sion fitness function satisfies the following gradient condition: 

d 2 f{x,y) d 2 f{x,y) 



dxidyj 

for all 



y=x 



dx j dy l 



y=x 



(19) 



Proof First of all, we note that in any case 



dsj(x) 
dx,; 



Therefore, 



d 2 f 



dx j dy l 



+ 



d 2 f 



y=x 



dyjdy, 



y=x 



D ^ x y = dsi{x) dsj(x) 



d 2 f 



dx q 



dxj 



dxjdy t 



d 2 f 



y=x 



dxidy. 



y=x 



(20) 



(21) 
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If the selection gradient is a gradient field, s(x) = VS(x), then 



(22) 



for all x, as claimed. Conversely, according to basic theory of differential forms 



(Spivak ( 1999 1), the condition = is equivalent to saying that the differential of 



the selection gradient is 0, i.e., that the selection gradient is a closed form. A closed 
form is exact if it is the differential of a function, and exactness of closed forms 



is captured by the first De Rham cohomology group (Spivak (1999)). For simply 
connected spaces this group is trivial, which implies that every closed form is exact, 
i.e., every closed form is the differential of a function, and hence a gradient field. 



Together, the two propositions above yield the following 

Corollary 2.3 Let f(x,y) be an invasion fitness function defined for a single species 
on a simply connected phenotype space of arbitrary dimension. Then the function 
f(x, y) is equivalent to the invasion fitness derived from a generalized Lotka-Volterra 
model for symmetric, frequency-dependent competition if and only if 



d 2 f(x,y) 



dxidyj 



d 2 f(x,y) 



y=x 



dxjdy % 



(23) 



y=x 



for all If this condition is satisfied, the generalized carrying capacity and com- 
petition kernel are given by expressions (17) and (18). 



3 Applications 



We illustrate the potential usefulness of the theory presented in the previous section 
with some examples. 
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3.1 Definition of frequency-dependent selection 

The first application concerns the conceptual issue of frequency-dependent selection, 
a fundamental and much debated topic in evolutionary theory. For adaptive dynamics 



models, it has been argued by Heino et al. ( 1999) that selection should be considered 



frequency-independent (or trivially frequency-dependent) if the resident phenotype 
x only enters the invasion fitness function through a scalar function describing the 
resident's effect on the environment (e.g., the resident population size). Otherwise, 
selection is frequency-dependent. If an adaptive dynamics model is given by an in- 
vasion fitness function whose selection gradient is a gradient field, this notion of 
frequency dependence can be made mathematically precise in the context of the gen- 
eralized competition model (15): selection is frequency-independent if and only if 
the generalized competition kernel a(x, y) given by (18) is a constant (equal to 1), 
i.e., if and only if 

/(*,y) = l-®, (24) 
K{y) 

where K(x) is the generalized carrying capacity (17). Because the generalized car- 
rying capacity is entirely defined in terms of the invasion fitness, this leads to a 
model-independent definition of frequency independence (and hence of frequency 
dependence). 



3.2 Universality of symmetric competition in 1 -dimensional phe- 
notype space 

The second application concerns the adaptive dynamics in 1 -dimensional phenotype 
spaces. In this case, any differentiable invasion fitness function trivially satisfies 
the gradient condition (19), i.e., the selection gradient s(x) is always an integrable 
function, and hence any 1 -dimensional adaptive dynamics model that is based on 
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a differentiable invasion fitness function is equivalent to the adaptive dynamics of 
symmetric logistic competition. We illustrate this by considering the invasion fitness 
function for a 1 -dimensional trait under asymmetric competition, 

a(x,y)K(x) 



f(x,y) = 1 - 



K(y) 



(25) 



where the competition kernel a(x, y) still has the property that a(x, x) — 1 for all x, 
but doc{x, y)/dy\ y=x ^ in general. Such models have been considered in the litera- 



ture (e.g. |Taper & Case] ( [T9^ , [KisdT| ( fT999) , pbebeli & Dieckmann| ( |2000l ), |Doebeh 
( 201 1] )) and reflect the assumption of an intrinsic advantage of one phenotypic di- 
rection, such as an intrinsic advantage to being higher or larger when competition 
between plant individuals is affected by access to sunlight. As in symmetric compe- 
tition models, the function K(x) is the carrying capacity function. 
In this case, it follows that 



s(x) 



9f{x,y) 



da(x,y) 



Let 



A(x) 



dy 



da{x,y) 



y=x 



dy 



+ 



y=x 



K'(x) 
K(x) 



(26) 



dy 



(27) 



y=x 



and consider the generalized carrying capacity 



- / A(x')dx' + \n(K(x)) 



K(x) = exp 



and the generalized competition kernel 

&{x,y) = (l-f(x,y)) 5-^-. 

K{x) 

Then K(x) > for all x, a(x, x) — 1 for all x, and da(x, y)/dy\ y=x 
Also 

a(x, y)K(x) 



(28) 



(29) 
for all x. 



f(x,y) = 1 - 



K(y) 



(30) 
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and hence the invasion fitness function for asymmetric competition is equivalent to 
the invasion fitness function of a symmetric competition model, in which the original 
asymmetry in the competition kernel is shifted onto the generalized carrying capacity 
function K(x). 

3.3 Cyclic adaptive dynamics in asymmetric single-species mod- 



Third, to further illustrate the implications of the gradient condition (19), we present 
an example of the adaptive dynamics that may occur when the gradient condition 
(19) is not satisfied. For this we consider an asymmetric competition model with a 
2-dimensional phenotype space defined by a carrying capacity function K(x) and a 
competition kernel a(x,y), where x = (xi,x 2 ) and y = (2/1,2/2) are 2-dimensional 
traits, as follows: 



We will not try to biologically justify this choice of functions, except to note that the 
competition kernel (32) has been chosen so that for i = 1,2, da(x, y)/dyi\ y=x 7^ 0. 
Biologically, this means that there are epistatic interactions between the trait compo- 
nents x 1 and x 2 , and as we will see, such interactions can be the source of compli- 
cated evolutionary dynamics even in a single evolving species. Note that the carrying 
capacity (31) has a maximum at (0, 0), that the competition kernel a given by (32) 
retains the property that a(x, x) — 1 for all x, and that for c\ = c 2 = 0, the com- 
petition kernel becomes symmetric. For c\ ^ or c 2 ^ 0, and in contrast to the 



els 




(3D 




(32) 



x exp (ci(xiy 2 - x 2 yi) + c 2 {x 1 (y 1 - x x ) + x 2 (y 2 - x 2 ))) . 
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1 -dimensional case for asymmetric competition kernels discussed above, the result- 
ing invasion fitness function f(x 1 , x 2 , y±, y 2 ) does not satisfy the gradient condition 
(19), and hence is not equivalent to the invasion fitness of a symmetric competition 
model. This can be seen by directly checking condition (19), or by considering the 
corresponding adaptive dynamics given by the selection gradients 



SilX 



S2(X 



9f{x,y) 



df(x,y) 



da(x,y) 



dy 2 



y=x 



y=x 



dyi 
da(x,y) 



+ 



dK(x) 



1 



dy 2 



dxi K(x) 
dK(x) 1 



+ 



y=x 



dx 2 K(x) 



(33) 
(34) 



It is easy to see that this adaptive dynamical system has a singular point at x* = (0, 0) 



(note that da(x, y) / dyi\ y=x=x * = for i = 1, 2), and that the Jacobian at the singular 
point x* has complex eigenvalues ±ici — c 2 . In particular, if c\ 7^ 0, the adaptive dy- 
namics has a cyclic component in the vicinity of the singular point x*, and it follows 
that if ci 7^ 0, the selection gradient s = (s 1 ,s 2 ) cannot be a gradient field because 
the gradient condition is not satisfied. In fact, the adaptive dynamics given by the 
selection gradients (33) and (34) can exhibit a stable limit cycle, as is illustrated in 
Figure 1 . Cyclic evolutionary dynamics are known to occur in the adaptive dynam- 



ics of multiple interacting species (Marrow et al. ( 1996), Law et al. ( 1997), Doebeli 



& Dieckmann (2000)), but the figure shows that such dynamics can also occur in 



the adaptive dynamics of a single species with multidimensional phenotype. The 
occurrence of cyclic evolutionary dynamics in a single species implies that the cor- 
responding adaptive dynamics model is not equivalent to a symmetric competition 
model (even if it is an asymmetric competition model, as in the example above). 
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4 Local universality of symmetric competition models 
for determining evolutionary stability 



So far we have considered the possible equivalence of general single- species adap- 
tive dynamics models in a given dimension to the adaptive dynamics of symmetric 
competition models. If the adaptive dynamics convergence to a singular point, the 
question of evolutionary stability of the singular point arises, and we may then ask 
whether the evolutionary stability of singular points of generic single-species adap- 
tive dynamics models can be understood in terms of the evolutionary stability of 
singular points in symmetric competition models. In this section, we show that in 
fact, the evolutionary stability of singular points of any generic single-species adap- 
tive dynamics models can always be understood by means of symmetric competition 
models. Consider an arbitrary invasion fitness function f(x, y) in the neighbourhood 
of a singular point x*, which we assume to be x* = without loss of generality. 
Because f(x*, x*) = as always, and 

df(x*,y)/dy\ y=x * = (df(x*,y)/dy 1 \ y=x *,...,df(x*,y)/dy rn \ y=x *) = by assump- 
tion of singularity of x*, the Taylor expansion of the function / in x and y in a 
neighbourhood of the point (x*,x*) has the following form: 



where g(x) is some function of x, A is a square matrix, and H is the symmetric 
Hessian matrix 



f(x, y) = g{x) + xAy T + -yHy T + h.o.t. 



(35) 




y=x=x* 



d 2 f(x,y) 
dyidy m 



y=x=x* 



\ 



H 



(36) 



d 2 f(x,y) 



l f{*,v) 

9y?n 



y=x=x* 



y=x=x 
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Here h.o.t. denotes terms of order > 2 in y, and x T and y T denote the transpose 
of the vectors x = (xi, x m ) and y = (yx,..., y m ), where m is the dimension of 
phenotype space. In generic models the evolutionary stability of the singular point is 
determined by the Hessian matrix H (Leimar (2009), Doebeli ( 201 1[ )). Here we call 



a model generic if the Hessian H is non-degenerate (i.e., has a trivial kernel). 
Now consider the modified invasion fitness function 

f(x,y) = g(x) + ^yHy T , (37) 

also defined in a neighbourhood of x*. Then the Jacobian J of the corresponding 
selection gradient s is simply the matrix H, which by construction is symmetric and 
hence has real eigenvalues. In particular, the linear dynamics defined by s has no 



cyclic component, and hence s is a gradient field (Spivak ( 1999)). It therefore fol- 
lows from Proposition 2.1 that the invasion fitness function / is equivalent to the 
invasion fitness of a symmetric competition model. On the other hand, the evolu- 
tionary stability of the singular point x* for the invasion fitness function / is also 
given by the matrix H, and hence is exactly the same as the evolutionary stability 
of the singular point x* for the original invasion fitness function /. This proves the 
following 

Proposition 4.1 Given any single-species invasion fitness function f with a non- 
degenerate Hessian matrix H, and a singular point x* of the corresponding adaptive 
dynamics model, then in a neighbourhood of x* , there is an invasion fitness function 
f derived from a symmetric competition model for which x* is a singular point, and 
such that the evolutionary stability of x* is the same for f and f. | 

It is important to note that while the evolutionary stability is the same for / and /, 
convergence stability of x* in the adaptive dynamics models defined by / and / (i.e., 
by s and s in the construction above) is generally not the same. We note that even if it 
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were, convergence stability generally depends on the mutational variance-covariance 
matrix (Leimar ( 2009[ )). Nevertheless, the proposition shows that the evolutionary 
stability of singular points of any single species adaptive dynamics model can be 
fully understood in terms of the evolutionary stability of singular points in symmet- 
ric competition models. We again note that evolutionary stability is not affected 
by the mutational variance-covariance matrix even in high-dimensional phenotype 
spaces (Leimar ( 2009| )), hence this result is important for generalizing results al- 
ready known for evolutionary stability in competition models, such as the finding in 
|Doebeli &I spolatov| ( |2010| ) that increasing the dimension of phenotype space gen- 
erally increases the likelihood of evolutionary branching in symmetric competition 
models. 



Conclusions 

We have shown that any model for the adaptive dynamics of a single species 
that is defined on a simply connected phenotype space and satisfies the gradient con- 
dition (19) is equivalent to an adaptive dynamics model for symmetric frequency- 
dependent competition. Specifically, expression (15) can be considered a "normal 
form" for any given invasion fitness function satisfying the gradient condition (19). 
The obvious advantage of having such a normal form is that results obtained for 
the normal form are general and hold for any adaptive dynamics model that can be 
transformed into this normal form. To illustrate this, we have shown that any single- 
species adaptive dynamics model in a 1 -dimensional phenotype space is equivalent 
to a symmetric competition model. In addition, we have shown that the evolution- 
ary stability of single-species adaptive dynamics models with arbitrary phenotypic 
dimension can be understood in terms of symmetric competition models. Attempts 
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at finding normal forms of invasion fitness functions have been made previously. For 



example, Durinx et al. (2008) showed that for every single-resident fitness function 
there exists a Lotka-Volterra competition model that has the same single-resident in- 
vasion fitness function. However, the interaction function was not partitioned into a 
frequency-dependent competition kernel and a frequency-independent carrying ca- 
pacity. The approach presented here appears to be at the same time simpler and more 
specific, which is probably due to the fact that purinx et al.\ ( |2008| ) dealt with the 
more complicated issue of normal forms describing the transition from monomor- 
phic to polymorphic populations. Our normal form (15) is simpler because it only 
considers selection gradients in monomorphic populations, and hence only requires 
the definition of the generalized carrying capacity function (17) and the generalized 
competition kernel (18). Even though our normal form is only valid for invasion 
fitness functions satisfying (19) and for the monomorphic resident population, for 
those conditions it is general because it holds globally, i.e., everywhere in phenotype 
space, rather than just in the neighbourhood of singular points. And it is more spe- 
cific because it disentangles the frequency-dependent and the frequency-independent 
components of selection. Generalized competition function without separation of 
frequency-dependent and frequency-independent components have also been con- 



sidered by |Meszena et a/.| ( |2005| ) for polymorphic populations. 

Essentially, the normal form (15) holds for any frequency-dependent adaptive 
dynamics model whose selection gradient is a gradient field, and hence whose dy- 
namics can be described as a hill-climbing process on a fixed landscape. Because 
such hill-climbing processes cannot capture oscillatory behaviour, this makes it clear 
that the normal form cannot hold for any adaptive dynamics model exhibiting cyclic 



dynamics, e.g. evolutionary arms races in predator-prey systems (Doebeli & Dieck- 



mann 



(2000), Doebeli (2011)). To illustrate this, we have given an example of cyclic 
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dynamics in a single- species adaptive dynamics model for asymmetric competition 
in 2-dimensional phenotype space, which therefore does not have the normal form 
(15). In general, the normal form can be applied to models of single species with 
high-dimensional phenotype spaces, and it is important to note that the generalized 
carrying capacity (17) and competition kernel (18) may be complicated functions 
in general. For example, the generalized carrying capacity may have multiple local 
maxima and minima (each representing a singular point of the adaptive dynamics), 
and the generalized competition kernel may have positive curvature at x = y (repre- 
senting positive frequency dependence). Accordingly, the adaptive dynamics result- 
ing from a normal form may exhibit repellers and dependence on initial conditions. 

The normal form (15) does not generally apply to the adaptive dynamics of mul- 
tiple species. For example, if a single species adaptive dynamics model has a normal 
form (15), and if that normal form predicts evolutionary branching, then evolutionary 
branching does indeed occur in the given adaptive dynamics model, but the normal 
form cannot be used to derive the adaptive dynamics after evolutionary diversification 
has occurred. This is because to derive the adaptive dynamics ensuing after evolu- 
tionary branching, one has to know the ecological attractor of coexisting phenotypes, 
i.e., one has to have explicit information about the ecological dynamics underlying 
the adaptive dynamics model, which the normal form does not contain. We leave it 
as a challenge for future research to derive normal forms for the adaptive dynamics 
of multiple interacting species. 
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Figure legend 



Figure 1: Example of cyclic adaptive dynamics in a single species with 2-dimensional 
phenotype space. The figure shows a numerical solution of the dynamical system 
dxi/dt = si(xi, x 2 ) and dx 2 /dt = s 2 (xi, x 2 ), where si and s 2 are given by (33) and 
(34). This system reflects the simplifying assumption that the mutational variance- 
covariance matrix (Leimar (2009), Doebeli (2011 1) is the identity matrix. Panel (a) 
shows the two traits x\ and x 2 as a function of time, and panel (b) shows the corre- 
sponding phase diagram, illustrating convergence to a limit cycle from two different 
initial conditions both inside and outside the limit cycle. Parameter values were 
o'k = Oa = 1, c\ = —1 and c 2 = —0.1. 
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